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Abstract 



'^ I We introduce a criterion how to price derivatives in incomplete markets, based on 

C ' the theory of growth optimal strategy in repeated multiplicative games. We present 

reasons why these growth-optimal strategies should be particularly relevant to the 
problem of pricing derivatives. Under the assumptions of no trading costs, and no 
restrictions on lending, we find an appropriate equivalent martingale measure that 
K^ ■ prices the underlying and the derivative security. We compare our result with other 

alternative pricing procedures in the literature, and discuss the limits of validity of 
^\j I the lognormal approximation. We also generalize the pricing method to a market 

^D ■ with correlated stocks. The expected estimation error of the optimal investment 

—J. . fraction is derived in a closed form, and its validity is check with a small-scale 

Q^ I empirical test. 
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1 Introduction 



The option pricing formula of Black-Scholes is arguably one of the most suc- 
cessful modern applications of quantitative mathematical analysis. Despite its 
great simplicity it gives prices of the most actively traded options on leading 
markets which are typically not off by more than a few percent. This has paved 
the way for an influx of advanced ideas from the theory of stochastic processes 
into finance. In the most general formulation of the no-arbitrage argument of 
Black-Scholes [11] and Cox et al [16] it was established by Harrison-Kreps [27] 
that a price system (under certain restrictions such as no trading costs) admits 
no arbitrage opportunities if and only if there exists an equivalent probability 
measure, with respect to which all discounted price processes of underlying 
and derivative securities are martingales. 

This theory is presented in several excellent monographs, see [17,29,32,21,28]. 
Nevertheless, these developments do not yield a fully general solution to the 
derivative pricing problem, since the equivalent martingale measure is not 
uniquely defined, except if the market is complete. If the market is incomplete 
there are many possible equivalent martingale measures, and each of these 
specifies a price system without arbitrage opportunities. 

In the lognormal model of Black-Scholes the market of one share and one 
bond is already complete. One can thus test its predictions by computing the 
volatility of the price process of the underlying security, plugging that into 
the Black-Scholes formula, and then comparing with observed market prices 
of a derivative: such tests were initiated already by Black-Scholes [10] and 
Black [12]. The volatility from historical data is obviously subject to mea- 
surement errors and statistical fluctuations. A more robust test is therefore if 
observed option prices can be fltted by a Black-Scholes formula with volatility 
treated as a free parameter. There are several recent reports that this is not 
so, see Rubinstein [36], and earlier investigations [31,35,15]. 

In addition to the lognormal model there exists a wider class of Ito processes 
in continuous time for which the market is also complete, and the derivative 
pricing problem can be solved only by no-arbitrage conditions. To distinguish 
which model in this class is the best approximation to observed data, and then 
to use it, is however already not a trivial task, see e.g [19,20,36,7,8]. 

For the rest of this paper we will consider the general case of an incomplete 
market, but we do keep the standard assumption of no market friction. We 
introduce a criterion to price derivatives in this situation. Stated differently, we 
advance an argument as to which out of many possible equivalent martingale 
measures is the appropriate one. We refer to this criterion as the Principle of 
No Almost Sure Arbitrage [1]: it is based on the theory of optimal gambling 



of Kelly [30]. We will consider an investor that gambles on investing in an 
underlying and a derivative security, and who looks for the growth-optimal 
investment strategies among such combinations. 

The paper is organized as follows: in section 2 we summarize and discuss the 
Kelly theory in a version appropriate for our purposes. We discuss also here 
the papers of Samuelson [37] and Merton-Samuelson [33]. In section 3, the 
technical core of the paper, we apply the Kelly theory to derivative pricing, 
and in section 4 we discuss the relations and differences to other proposed 
alternative derivative pricing procedures. 

It is straight-forward to see that in complete markets our procedure agrees with 
the standard ones, e.g. with the models of Cox-Ross-Rubinstein and Black- 
Scholes. In section 5 we discuss the lognormal limit of our procedure, the 
limits of its applicability, and that it is not exactly equivalent to the Black- 
Scholes formula. In section 6 we study the expected scatter of the optimal 
investment fraction and report a small-scale empirical test of pricing options 
on the Swedish OMX index. In section 7 we generalize the method to a market 
with correlated stocks. In section 8 we summarize and discuss our results. 

An earlier version of this paper was circulated in 1998 [2]. A longer sequel 
was also made available in the electronic archive [1], lacking however in the 
version accepted for publication, material corresponding to sections 6 and 7 
below. The presentation in the present paper is significally different from [1] 
and contains a fuller discussion of the implementation of the model and the 
economic background. 



2 The Kelly problem of Optimal Gambling 



The theory to be exposed in this section is due to Kelly [30], who was looking 
for an interpretation of Information Theory of Shannon [43] outside of the 
context of communication, and to the treatment of the St. Petersburg Paradox 
by Bernoulli [9]. For a recent review of growth-optimal investment strategies, 
see Hakanson-Ziemba [25]. 

Consider a price movement of stock or some other security which is described 
by 

S',:+i = UiSi (1) 



where time is discrete. Si is the price at time i and the -Uj's are independent, 
identically distributed random variables. The assumptions on the Mj's can be 
relaxed. 



Consider now an investor which starts at time with a wealth Wq, and who 
decides to gamble on this stock repeatedly. Suppose that the investor chooses 
to commit at each time a fraction Ig of his capital in stock, and the rest in 
a risk-less security. We set the risk-less rate to zero. A non-zero risk-less rate 
corresponds to a discount factor in the definition of the share prices, and can 
be accounted for by a redefinition of the Mj's. The investor will hence hold at 
time i a number IgWi/Si of shares, and his wealth at successive instants of 
time follows a multiplicative random process 

W^+i = (1 - ls)W^ + huiWi = (1 + h{u, - 1)) Wi (2) 



In the large time limit we have by the law of large numbers that the exponential 
growth rate of the wealth tends with probability one to a constant. That is, 

A(y = hm i log ^ = E^[log(l + /.(«-!))]. (3) 



The optimal gambling strategy in the sense of Kelly consists in maximizing 
\{ls) in (3) by varying Ig. The solution must be unique because the logarithm 
is a concave function of its argument. Which values of Is are reasonable in our 
problem? First, the optimum /« must be such that l + ls{u — 1) is positive on 
the support of u. Second, we must decide if the investor is allowed to borrow 
money. In the original formulation of Kelly he is not, but here it is useful to 
assume that the investor may do so at a risk-less rate. Furthermore we also 
assume that the investor can go short i.e., it is possible to borrow stocks. 
Hence we allow Is to take any finite positive or negative value, and look for 
the maximum of \{ls)- 

The desired strategy is hence the only finite Is which solves 

—r — i,=i,* - t. , - U (4j 



and the maximum realized growth rate is 

\* = E^[\og{l + ls*{u-l))] (5) 



Let us look more closely at (4) if u can take the values Vn with probabilities 
Pn- Let us write 

Pn /„N 



such that (4) reads 

^g„(t;„-l) = (7) 



From 1 + ls*{vn — 1) being positive it follows that the g„'s are all positive, and 
from (7) it follows that the sum of the g„'s is one. The g„'s thus form a new 
set of probabilities, and, again by (7), with respect to these the stock price 
process is a martingale. 

It is worthwhile to discuss this point at some length in light of the well- 
known papers of Samuelson [37] and Merton-Samuelson [33]. These authors 
demonstrated that if the risk/return preferences of an agent is described by 
an expected utility, then the growth-optimal strategy does not maximize that 
utility, except if the utility is logarithmic, and the price process is similar 
to (1). 

It is a mathematical fact that the growth rate (3) holds with probability one. 
If a gambler would choose any other strategy in Ig but that which maximizes 
A then he would in the large time limit almost surely end up with an expo- 
nentially smaller capital. 

In the counter-examples of Samuelson-Merton the dominant contribution to 
non-logarithmic utility comes from events with asymptotically (in time) van- 
ishing probability. When sets of measure zero are involved, the relation be- 
tween sample and ensemble averages becomes dehcate. Suppose some agents 
want to maximize a non-logarithmic utility, and that they choose strategies ac- 
cordingly. If they do so, and we compare with them using the growth-optimal 
strategy, they would almost surely end up with less utility according to their 
own criterion. Only if there would be an ensemble of infinitely many agents, 
all using this same strategy, then some few would actually end up with higher 
utility. And, in fact, so much higher that the averaged utility over the ensemble 
would be increased. This is the sense in which other strategies can be 'better' 
than the growth-optimal strategies. Therefore, Merton-Samuelson give prefer- 
ence to ensemble averages over sample average. This is a problematic point, 
which has apparently not received much attention in the economic literature. 

A simple counter-argument against the growth-optimal criterion is that time 
in economical problems, in the sense considered here, is perhaps often not 
very long. However, this is essentially a quantitive question, which has to be 
decided case-by-case. For recent discussions of characteristic times in the Kelly 
problem with transaction costs, see [3,6,42] 



3 Fixing the price of the derivative 



Let us now consider the stock of section 2 and one derivative security on the 
same stock. The argument generalizes at every step in an obvious way to an 
arbitrary number of derivatives, but to keep the notation simple we will write 
out the formulae for one derivative. To begin, let us write the present share 
price at the time i Si, and let us assume that the realized value of the derivate 
at the next future instant of time is Cj+i (5*4+1) = Ci+i{SiUi). The unknown 
parameter is the present price of derivative, which we will write Cj. We will 
further suppose that this situation is repeated over and over again. That is, we 
assume that the investor finds himself many times in the situation that if he 
keeps some money Wi in stock it will be worth UiWi, while if he keeps the same 
money in the derivative it will be worth Ci+i{SiUi)Wi/Ci{Si), with the same 
random variable Ui entering into both expressions. The simplest example is 
one time step before expiration time of a call option with a present share price 
St-1 and a strike price K; in this case one has Ct{St) = max{ST — K,Q). 

If at the i'th stage in this gamble the investor commits a fraction Ig of his 
wealth in shares and a fraction Id in the derivative, this means IsWi/ Si shares 
and IdWi/Ci derivatives. After the random variable Ui has taken a value, the 
wealth is changed to 

lV,+i =(l + h{u, - 1) + u^ -i)]w. 



a 



If the investor plays this game many times, his wealth will almost surely grow 
at an exponential rate (over these games) which is 

A(/„ la; Q) = ^^[log (l + /,(«-!) + /,(^ - 1))] (9) 



Let us now assume that we compute (9) and that we find an optimal growth 
rate 

\*{a) = Max,,,;,EP[log (1 + h{u - 1) + /,(^ - 1))] (10) 



and that this is realized at fractions /* and l^- If /^ is larger than zero, then all 
operators would like to buy the derivative. Hence its present price Ci would 
tend to rise, and its rate of return would fall. The fraction l^ would thus tend 
to fall. If, on the other hand, /^ would be less than zero, then all operators 
would want to go short on the derivative, its present price would fall, its rate 
of return would rise and l^ would tend to rise. The only value of C, in which 
the market can be in equilibrium is therefore the one such that /^ is zero. This 



statement was referred to as the Principle of No Almost Sure Arbitrage in 
[1]. We remark that from the technical point of view, this procedure, where 
l^ is taken to vanish, is closely similar to the Samuelson's [38] methodology of 
warrant pricing in the 'incipient case', see Merton [32, chapter 7]. 

Ci and the fractions Ig and l^^ thus satisfy two equations. One is 

d\{ls,ld;Ci) M-l 



which is identical to (4) of the Kelly model, and therefore determines /*. The 
other is 

^ rfA(/„/,;a), ^,, (%f-l) 



\ls=ll,li = Q,C^=C' 



If we now compare with (6) and (7) we see that (12) can be written 

C* = Y.(lnCr+l{S,Vn) (13) 

n 

which states that the derivative price process is a martingale with respect to 
9- 

We have constructed the equivalent martingale measure q over one elementary 
time-step. By compounding we can construct, from T convolutions of q dis- 
tribution, the price distributions over many time-steps Qt{.St\Sq) . The price 
of, say, a European call option with expiration time T at time t, can thus be 
written in the standard manner as 

C(i, Si) = E^^-' [C{T, St)] C{T, St) = max{ST - K, 0) (14) 

where the expectation value is taken with respect to the many time-steps 
distribution QT-i{ST\Si) over share prices St at time T, conditioned by the 
share price having been Si at time i. We have taken the risk-less rate to be 
zero. A non-zero risk-less rate can be put back in as a discount factor in the 
share prices, which leads to a different definition of the -u's, the returns over 
elementary periods of time, and hence to different g's. 



4 Comparison with other derivative pricing prescriptions 



The pricing procedure proposed in section 3 is based on a particular choice of 
an equivalent martingale measure. It therefore follows that this procedure is 



arbitrage- free. It also follows that it must agree with no-arbitrage pricing in 
complete markets. 

It may nevertheless be instructive to carry through the calculations of sec- 
tion 3 explicitly for the Cox-Ross-Rubinstein model. We hence assume that 
the discounted share price can go up by a fraction u/r with probability p, 
down by d/r with probability (1 — p), and that u > r > d, where r is the 
risk-free interest rate. Equation 11 is then solved for /* as 

»_ rp r{l-p) p _ r-d 

' r-d u-r '^ l + /*(f-l) u-d ^^ 



and this reproduces the Cox-Ross-Rubinstein measure g, independent of p. 
The maximum growth rate however still depends on p, and is 

A* = — plog (1— p)log (16) 

q l-q 



The Black-Scholes model can now be derived as the continuous-time limit of 
the Cox-Ross-Rubinstein model. 

The growth rate (16) is the Kullback contrast of the measure q with respect 
to the measure p in the dichotomic case. It is easy to see that this is the 
general form of the solution, and that the measure q can be obtained as follows: 
Consider a security with price process determined by a measure p, and consider 
all measures q with respect to which the price process of the security is a 
martingale. Then the particular measure q which we compute from (11) (see 
also equations (6) and (7)) is the one that minimizes the Kullback contrast 

dM = E^[-\og^] = -J2p^\og^ (17) 

Q n Qn 



Our proposal therefore coincides, over one time step, with one of the two 
proposals recently put forward by Stutzer [44]. Stutzer's formulae differ from 
ours in that we posit the minimization of (17) under the martingale constraint 
for the elementary process over each time-step, while Stutzer applies (17) 
directly on the probability distribution over a finite time interval. A similar 
approach minimizing the relative entropy distance to calibrate a pricing model 
has been suggested by Avellaneda et al [4]. 

It is also interesting to remark that if we perform the quadratic approximation 
of the logarithm in (3) around the expected return on capital, EP[Wi^i/Wi], 



and then look for the optimal portfolio, we obtain 



q{u) = 1^1 - ^^^^{u - 1 - fi)j p{u) ; fi = E^[u-l] {li 



Equation (18) reproduces the least-square option pricing procedure of Follmer- 
Sondermann [24], Schal[39] and Schweizer [40]. Another approach to option 
pricing based on global least-square minimization was recently proposed by 
Bouchaud-Sornette [14], see also [13], but has now been proven by Wolczyn- 
ska [46] and Hammarlid [26] to give the same least-square solution, see also 
Schweizer [41]. The empirical tests of [14] and [13], which pertain mainly to 
the case of small n, hence also validate our theory. 

A conceptual problem of least-square option pricing, first shown by Dybvig- 
Ingersoll [22], is that it can assign negative prices of some state-contingent 
claims. This may lead to negative prices of derivatives with strictly non- 
negative pay-off, and hence to arbitrage opportunities. In (18) this happens 
when u is sufficiently large, provided /i is positive. From our point of view 
this result is clear. Negative prices appear only at capital values for which 
the quadratic approximation is a decreasing function. It is precisely because 
quadratic utility is not monotonically increasing that its maximization only 
gives a pseudo-martingale measure, which can take both positive and negative 
values. Our theory is practically equivalent to the least-square method when 
that one has no problems, i.e., when we look at small deviations from the most 
probable realizations of the capital returns and when /i is small. In the oppo- 
site limits our theory deviates from least-square minimization, always prices 
derivatives by a proper equivalent martingale measure, and the problem of 
possible negative prices does not appear. 



5 On the lognormal limit 



The distribution Qt{St\So) is constructed by compounding the measure q of 
one time-step. If u only takes only a finite number of values, then Qt{St\So) 
can be written as a multinomial expansion. However, if the number of time 
steps is large, then the multinomial expansion is unwieldy, and some other 
approach must be found for analytic and numerical work. 

The asymptotic price distribution is of course determined by the behavior of 
aggregated logarithmic returns, and this is intimately linked to the lognormal 
distribution, its applicability and its limits. Let us assume that E'^[(\ogu) ] 
is finite. Then, for the most probable realizations, Qt{St\So) is well approxi- 



mated by a lognormal distribution, i.e. 

^ ^c \c\ niiL^)ic \c\ ^ {\og{ST / Sq) — XT) 



where 

\ = E''[\ogu] ■ A^ = E\{\oguf] - {E''[\ogu]f (20) 

For large T only low-probability events are not distributed according to (19). 
For the treatment of such large, but very rare, events the proper mathematical 
setting is the theory of large deviations, see Varadhan [45]. In [1] we discuss 
the use of large deviations theory to price derivatives related to rare events. 

We note that Qt{St\Sq) is by construction a martingale measure, while, for 
the log- normal approximation, 

EQ^""^[5^] = |Q^^^)(5T|5o)5Trf^T = ^oexp(A + ^A2)T . (21) 

In general A + |A^ ^ 0. That is, in the lognormal approximation one looses 
the martingale property. 

Our lognormal approximation does not agree with a straight-forward appli- 
cation of the Black-Scholes formula. We recall that the Black-Scholes theory 
also involves lognormal distribution. In the parameterization of (19) it is char- 
acterized by drift coefficient A^^ and volatility A^^g, related through 

Xbs = -^^Is , A|^ = E^iloguf] - {E^[\ogu]r (22) 



It is easy to understand that if u has support concentrated around one, then 
the difference between Qj< and Qrp is small. In the continuous limit the 
differences between Qt, Qt and Qrp disappear altogether. 



On the other hand, the risk-neutral price distribution can actually be close to 
lognormal but still significantly different from Black-Scholes. In an extreme 
situation of very large l^ar^[-u], Qj^ is only close to Qt in a very limited 
region. As an example we show in figure 1 the case where p{u) is a truncated 
Levy distribution. The distribution Qt has been computed by a Monte Carlo 
procedure. Here it is obvious that one has to consider the corrections to the 
lognormal to get a good approximation to Qt- 

We observe that if u has support around one the predictions of the Black- 
Scholes theory and the lognormal approximation are both rather close to the 
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Black-Scholes lognormal distribution - 

Lognormai approximation - 

Monte Carlo simuiation - 




Fig. 1. The probability distribution is a discrete approximation to a Levy distribu- 
tion. The relative price increments can take the values Un = uoa'"", with uq equal to 
0.4, a equal to 1.2 and the index n ranging from zero to N, which in this example 
has been set to 30. The probability of an elementary event, p{un), is CnUh , with 
P equal to 0.5 and Cm a normalization constant. The Monte Carlo simulation was 
performed with 20 million trials. 

exact formula: this is the case in figure 2. In figure 2 we show the price of 
a European call option as a function of strike price where the probability 
distributions are: the increment u takes three discrete values 0.8, 1 and 1.2 
with probabilities 0.2, 0.3 and 0.5. 

On the other hand, in general, both Black-Scholes formula and the lognor- 
mal approximation can be far from the exact formula. Large deviations have 
therefore to be taken into account to compute properly the expectation value, 

see [1]. 



6 The experimental error on the optimal /* 



The realistic application of the theory we have presented depends on the pos- 
sibility of estimating the distribution p{u) from real data. In fact, p{u) enters 
into the construction of the q{u) both directly and with the optimal /*, which 
is a functional oi p{u). Therefore, it is very important to estimate I* together 
with the magnitude of its error. 

The information which is possible to obtain from the market is a series of stock 
prices. So, ■ ■ ■ Sh, where H is the number of data. We consider the time interval 
between observations to be large enough to ensure the time independence of 
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Fig. 2. Option prices C as function of the strike prices K for T = 30. The relative 
price increments over one elementary step in the process are, that u takes three 
discrete values 0.8, 1 and 1.2 with probabilities 0.2, 0.3 and 0.5. 

the returns Ui = Si/Si^i. 

In order to estimate experimentally the true distribution p{u) we consider a 
partition of the support of the distribution in M intervals of width A, the 
requirement being that M < H. Then we estimate the probability in one of 
the interval by the ratio between the number of observed data in the interval 
and the total number of data H. In this way one creates a histogram of the 
approximating probability Pa,h{u), the deviation of which from the 'true' 
probability p{u) decreases as H becomes larger and A smaller. Let us call this 
deviation 



SPa,h{u) = PaAu) - P\ 



u\ 



(23) 



Because an observed data is either in an interval or not, the number of ob- 
served data in an interval follows a binomial distribution. This observation 
immediately leads to 



\/ar(5PA,//H) = Pa,//(w)(1 - Pa,//HA) 



(24) 
(25) 



To estimate the true /* with l*^ using the approximating probability Pa,h{u) 
we have to solve the equation 



Pa,//(^)(^-1) 
l + lhiu-l) 



du = . 



(26) 
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The error on the optimal fraction invested in the stock is the difference: 

6l = l*H-l*, (27) 

where /* is the true optimal fraction obtained from the distribution p{u). 

A Taylor expansion up to the second order around /* of (26), and the above 
considerations about the error in estimating the probability give after some 
calculations 



where 

\J (l + /*(«-l))2 i ^ ^ 



where it has been used the fact that for a sensible partitioning of the sample 
space PA,_f/(w)A ^ 1. 

The error on /* can also be estimated from the historical distribution using 
bootstrap technique [23], i.e. drawing with replacement a set of J new samples 
of the same size H of the original one from the empirical distribution, and 
computing a new /* for each bootstrap sample. In figure 3 we compare the 
variance of this distribution with the quadratic estimation (28). The agreement 
is very good for different values of the (T„ = JVar[u\. 

According to equation (28), it is evident that in most cases the error of /* is 
large. A non-intuitive result is that it is easier to determine /* with accuracy 
for distribution with larger variance because of the a in the denominator. This 
fact can be understood as follows. A stock with small fluctuations allows for 
a bigger fraction be held against the interest rate to make money, but just a 
small change in the distribution will make a big change in the optimal fraction. 
On the contrary, for a wildly fluctuating stock the fraction invested is smaller, 
and a little change in the distribution will not have much effect on the optimal 
portfolio. 



6.1 Data analysis 



Empirical study is needed to verify the growth optimal pricing procedure here 
discussed. To verify the model we looked at the European call option on the 
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Fig. 3. Comparison between the error on I* estimated by equation (28) and by 
the bootstrap technique. We plot l/{Haf) vs o"^, where the ai are obtained by the 
bootstrap technique with J = 200 samples and the distribution of m is a gaussian 
with E[u] = 1 and standard deviation o"„. 

OMX-index, traded at the OM exchange in Stockholm, during the period 2 
January 1995 until 22 September 1995. We consider closing prices of each day. 
The time period between trading is one day. The time series contained 183 
observations, providing a sample of Ui of size (data points) 182. The empirical 
distribution of Ui of the OMX-index is built by the fraction of the price of the 
OMX-index today and tomorrow, Ui = S'j+i/S'j. Weekends and other holidays 
are treated as non-existing. 

The thirty day interest rate was as used discount factor r to compute the 
discounted daily returns, {ui/ri}}^^. During the period the interest rate fluc- 
tuate between, 7.15 — 9.83% of yearly yield. The sample space was divided 
into M = H cells and the optimal /* determined using the Newton-Raphson 
algorithm [34]. Because of the concave of the logarithm the convergence is 
very fast, the optimal /* for the OMX-index is equal to 14.46. 

The distribution of /* is obtained using bootstrap technique with a set of 200 
new samples. The standard deviation of /* is 8.6 and can be compared to 
8.7 obtained estimating the standard deviation using formula (28). From the 
empirical distribution p the empirical g-distribution can be created and from 
g-distribution the empirical Q-distribution can be derived. 

The empirical g-distribution is used to estimate parameters of the lognormal 
approximation (20) A = 1.9 ■ 10~^ and A^ = 7.3 ■ 10~^. We determine, using 
the bootstrap technique, the distribution and standard deviation for A and 
A^, ax = 3.6 ■ 10-6 and aA^ = 3.9 ■ 10"^ respectively. 
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Fig. 4. Experimental option prices against theoretical ones for the call options of 
OMX with maturity date less or equal then 10 working days in the period between 



the 2 of January and the 22 of September 1995. The theoretical prices are 
generated with the probability QTiSr/So) using the Monte Carlo discussed in this 
section with 10^ trials. The linear fit with ax + b gives a = 1.003 it 0.004 and 
b = O.iblO^^. As discussed in the text the major source of errors is on the estimation 
of /* due to the low data frequency considered. 

These parameters correspond, inside the experimental errors, to the Black- 
Scholes model. The consequence is that in this case lognormal approximation 
of the Kelly pricing scheme is practically equivalent to Black-Scholes model. 

The comparison between empirical data and theoretical one is shown in fig- 
ure 4, for the call options with maturity date less or equal then ten work- 
ing days (two weeks). We observe a coincidence inside the errors. Finally we 
observe that to be able to tell a difference between classical way of pricing 
derivatives and growth-optimal pricing procedure either the price movement 
of the underlying has to have 'wilder' fluctuations, or one has to have access 
to more high-frequency data. 



7 The general case with correlated stocks 



In order to study the effect of correlated stocks, the market in this section is 
assumed to have L stocks {S''}j^^^ and D^ derivatives {C^^'''\S!^)}j^i written 
on them. To get transparency of the principle we discuss the case of just two 
stocks L = 2 with one derivative each and a risk free interest rate equal to 1. 
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The conclusion of this section will be given in a general notation. 
The features of the stock price movements are given by: 

qik) 

-f^ = Sr ^ = 1^2. (30) 



The random variable u\ is assumed to be correlated on other stocks. 

An investor of thgs market commits a fraction of his capital l^ in stock k and 
a fraction dj^k in the derivative j. The exponential growth rate of the capital 
is: 



A(/i, I2, di^i, di^2) 



(1.2) 

1 



.2) 



log (1 + /i(tx(l) - 1) + h{u^^^ - 1) + rfl,l (^ - 1) + rfi,2 (1^ 

According to the general equlibrium argument, one obtains the two equations: 

E« [(«('=)- 1)] = A; = 1,2 (32) 

and 

cF\sf^) = i^ncftf (5ffi)] J = 1, 2 (33) 

where we have defined, as in section 3, a new probability measure q: 



Qh 



(2) _ -^ Pfv 



tr 1 + iiiuc^'^^ - 1) + /^(m(™'2) _ 1) ^"^ ^ 



where the index h and m stands for some realization of the random variables 
u^^^ and -u*-^-* with the joint probability phm- 

The solution of equations (32) gives the optimal fraction II to invest in each of 
the different stocks S'' and equations (33) give the option prices as the expected 
value under the measure (34). The martingale measure may be compounded 
so that the option is priced by using the pay out function at expire date: 

^p-,fc)(^W) ^ eQ[4^»(4^))] (35) 

where the Q denotes the compounded martingale measure of q. 
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(31) 



The case of independent stock price movements is of particular interest from 
a theoretical point of view. The price of the derivative on an underlying stock 
and even the hedge position do not depend on the other stocks in the 'classical 
models', i.e. complete markets and risk minimization pricing procedure. This 
result is a consequence of the property 

Qh^ = (lh (36) 



(2) 

where qu and q\ are the martingale measures of one and two shares respec- 
tively. 

In the case of a complete market it is a known fact that the martingale measure 
Qh is unique [18], then it can not be anything else but the measure (6) found 
by optimizing the growth rate for just that particular stock 

^'^ = i+/*(r(M)-i) • (37) 

The property (36) is also preserved in the risk minimization pricing procedure 
but not in the general case of an incomplete market. In this situation the 
associated martingale measure (34) depend on the number of stocks traded in 
the market 

(2) _ Y^ TPjJ}: -L l'XR\ 

lh -Ph-2^^^ ;*^^(,,i) _ ^) ^ ^*(^(™,2) _ 1) ^ ^^ • ^"^^^ 



In figure 5 we show this difference on the option prices for the trichotomic 
case. 

This fact is clear in an incomplete market with large number L of stocks and 
the returns of these stocks are independent identically distributed random pro- 
cesses. The optimal fraction invested in each stock is 0{1/L) and the excess 
rate of return becomes almost sure. To avoid arbitrage opportunities the ex- 
pected excess rate of return n must vanish and the measure pricing the option 
becomes close to the observed probability i.e, C ~ E^[Ct]- Let us notice that 
this corresponds to a generalization to a non gaussian p{u) of the Bachelier 
theory [5]. 



8 Conclusion 



We have in this paper introduced a new criterion to price derivatives in in- 
complete markets derived from the theory of optimal gambling strategies in 
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Fig. 5. Option prices C rescaled by the initial value of the stock 5*0 vs. \og{K / Sq) 
for T = 30 and the returns of figure 1. We compare the exact prices with one and 
two asset in the market with the Black-Scholes prescription. 

repeated games. The criterion say that one should not be able to construct 
a portfolio using derivatives that grows almost surely at a faster exponential 
rate than using only the underlying security. 

A corner-stone in modern derivative pricing theory is that discounted prices 
(in friction-less markets) can be expressed as expected pay-off, expectations 
taken with respect to equivalent martingale measures. We use the criterion to 
decide which out of many possible measures to choose. We therefore do not 
need to assume a complete market to fix unambiguously a price. 

In section 4 we compared our approach to other procedures proposed in the lit- 
erature, and in section 5 we discussed the lognormal limit. In our construction 
we determined first the equivalent martingale measure over one elementary 
time-step. By composition we can then get the risk-neutral distribution at 
finite times. In a large class of models this compounded probability tends to 
the lognormal form, at the points close to the maximum, where probability is 
high. One could then conclude that we just rederive the Black-Scholes model, 
except for rare events. This is however not true, because, for us, the condition 
that the risk-neutral distribution is a martingale involves also the distribution 
of rare and large events. These are not lognormal. In fact, the best lognormal 
approximation does not by itself have to be a martingale, and in general differs 
from Black-Scholes. 
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